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Here we have systematically studied the cooperative binding of substrate molecules on 
the active sites of a single oligomeric enzyme in a chemiostatic condition. The average 
number of bound substrate and the net velocity of the enzyme catalyzed reaction are 
studied by the formulation of stochastic master equation for the cooperative binding 
classified here as spatial and temporal. We have estimated the entropy production for 
the cooperative binding schemes based on single trajectory analysis using a kinetic 
Monte Carlo technique. It is found that the total as well as the medium entropy 
production show the same generic diagnostic signature for detecting the cooperativity, 
usually characterized in terms of the net velocity of the reaction. This feature is also 
found to be valid for the total entropy production rate at the nonequilibrium steady 



state. We have introduced an index of cooperativity, C, defined in terms of the 
ratio of the surprisals or equivalently, the stochastic system entropy associated with 
the fully bound state of the cooperative and non-cooperative cases. The criteria of 
cooperativity in terms of C is compared with that of the Hill coefficient and gives 
a microscopic insight on the cooperative binding of substrate on a single oligomeric 
enzyme which is usually characterized by macroscopic reaction rate. 
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I. INTRODUCTION 



Conventional thermodynamics at or near equilibrium needs serious modification to accom- 
modate the events of single molecular processes as well as nano-systems which are generally 
in states far away from equilibrium^^. The single molecule study is very important in bio- 
logical systems because most of the processes in cell are taking place on the level of a single 
or a few molecules. The non-equilibrium feature is mainly developed within a cell due to 
the mechanical or chemical stimuli which runs the metabolism through the driven chemical 
reactions^^ 2 -. Quantitative measure of fluctuations^ 3 - - — in small systems are possible over 
short periods of time that allow the thermodynamic concepts to be applied to such finite 
systems. A crucial concept in the statistical description of a nonequilibrium small system is 
that of a single trajectory or path^^ 3 - 1 ^ 1 ^ and defining the entropy of the system for such a 
single trajectory allows one to formulate the second law of thermodynamics at the stochas- 
tic level 2 ^ - — . The trajectory-based entropy production 2 ^ - — has been successfully applied to 
various systems, for example, single bio-molecular reactions^, chemical reaction networks 22 
and driven colloidal particles^ 3 -. 

Enzyme kinetics is a very important process in cellular metabolism where nonequilibrium 
feature is developed due to the imbalanced chemical reactions and the presence of chemio- 
static condition prevents the reaction system to attain equilibrium 2 ^ - —. In a chemiostatic 
condition, substrate and product are maintained at constant concentrations by continuous 
influx of the substrate and withdrawing the product from the system. Under this condition, 
the reaction system reaches a nonequilibrium steady state (NESS) 25-27 which is character- 
ized by a non-zero total entropy production rate. Single molecule enzyme kinetics 2 ^ 3 ^ is 
theoretically studied using the stochastic master equation approach^ 2 - as well as by the 
stochastic single trajectory analysis^^ 2 -. Now most of the enzymes found in enzymology 
are oligomeric in nature consisting of two or more subunits usually linked to each other 
by non-covalent interactions 33 . Possibility of interaction between the subunits during the 
substrate binding process can give rise to different cooperative phenomena^ 3 - - —. Positive 
cooperativity is said to occur when the binding of one substrate molecule with a sub- 
unit increases the affinity of further attachment of the substrate to another subunit^ 3 - 1 ^ 1 ^ 7 -. 
In the case of negative cooperativity, attachment of a substrate molecule to one sub- 
unit decreases the tendency of further attachment of the substrate molecules to other 
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subunits^ME. These types of cooperativity based on the affinity of the substrate bind- 
ing belong to the class of allosteric cooperativity— ~—. There is another type of coopera- 
tivity, termed as temporal cooperativity 4 ^, reflected in the zero-order ultra sensitivity of 
the phosphorylation-dephosphorylation cycle which is shown to be mathematically equiva- 
lent to the allosteric cooperativity^ 1 -. Beside allosterism, cooperativity has been studied in 
monomeric enzymes with only a single substrate binding site. This has led to the important 
concept of hysteretio^ 1 ^ 1 ^ and mnemonic enzymes 4 ^ 4 ^. These two types of enzymes show 
the cooperativity phenomena due to the slow conformational disorder of the active site^. 

In this paper, we have studied the entropy production in the kinetics of a single oligomeric 
enzyme which shows cooperativity with respect to the substrate binding. Here we have clas- 
sified the cooperativity phenomena according to the nature of the different substrate binding 
mechanisms, namely, sequential and independent, as detailed by Weiss 44 . In sequential bind- 
ing, the adjacent sites of the oligomeric enzyme are successively occupied by the substrate 
molecules. So the substrate-bound states of the system are actually adjacent in space and 
hence we denote the cooperativity arising out of this binding protocol as the spatial co- 
operativity. For the sequential mechanism, the first binding site, i.e, the first subunit of 
the oligomeric enzyme must be filled in order for the second site to become occupied by 
the substrate, as if the substrate molecules have been stacked on top of each other at their 
binding sites 44 . This type of binding can be relevant to an ion transporter, such as the Na-K 
pump 45 . The other class is called temporal cooperativity which can occur due to the inde- 
pendent binding of the substrate molecules to any one of the subunits at a particular time 
without any specific spatial arrangement. Here the substrate-bound sites are not physically 
neighboring in the enzyme 4 ^ but the global state of the system is defined in terms of the 
total occupancy of the overall sites at a particular instant of time. This type of binding 
can be observed in multimeric proteins with individual binding sites located on different 
subunits, such as ligand gated ion channels or ligand gated enzymes^ 4 -. Here we have theo- 
retically studied the cooperative behavior solely from the viewpoint of the substrate binding 
mechanism and not in terms of the active and inactive enzyme conformations or the actual 
structural details of the enzyme that can lead to such mechanisms^ 4 ^. To study the bulk 
kinetics of allosteric enzymes Monod, Wyman and Changeux (MWC) in 1965 and Koshland, 
Nemethy and Filmer (KNF) in 1966 put forward models to account for cooperative binding. 
Generally 'Sequential' is used as a term for a classical distinction between multi-step bind- 
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ing models, as for instance differentiating between the KNF and MWC models in terms of 
the variation of the substrate binding rates in each successive step. So the term sequential 
in KNF model should not be confused with the term sequential used in our approach^^. 
We have constructed the master equations for each class of substrate binding. Time evo- 
lution of such cooperative systems can be described by suitably applying a kinetic Monte 
Carlo technique^ 1 ^. Here we have applied this algorithm to calculate the total, medium 
and system entropy production along a single trajectory for such cooperative systems as a 
function of the substrate concentration over a time interval where finally the system reaches 
a nonequilibrium steady state (NESS) and then determined the ensemble average quantities 
over many realizations of such trajectories. We show the correspondence between the evolu- 
tion of the total and the medium entropy production with the average substrate binding and 
net velocity of the reaction in the context of detection of the cooperative behavior. Similarly 
this correspondence is also studied for the total entropy production rate at the NESS. The 
system entropy production is thoroughly studied in terms of the substrate binding proba- 
bilities for the different classes of cooperative systems considered. We have introduced a 
cooperativity index, C defined in terms of the stochastic system entropy to understand the 
nature of the cooperativity. 

Layout of the paper is as follows. In Section II, we have given the master equations 
and their steady state solutions to describe the spatial and temporal cooperative binding 
mechanisms and the corresponding entropy production rates. In Section III, numerical 
results of entropy production and cooperative kinetics is discussed. In Section IV, we have 
discussed on measures of cooperativity and introduced an index of cooperativity. Then the 
paper is concluded in the Section V. 

II. COOPERATIVE BINDING, MASTER EQUATION AND ENTROPY 
PRODUCTION RATE 

In this section, we have first classified the cooperativity of a single oligomeric enzyme 
on the basis of the nature of the enzyme-substrate binding and then proposed a stochas- 
tic description for each class in terms of a one-dimensional random walk problem. Here 
we have provided a master equation approach for the description of spatial and temporal 
cooperativity which is suitable for the calculation of entropy production. 
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A. Classes of cooperativity: spatial and temporal substrate binding 



Here we have considered that the substrate molecules can bind to the subunits of the 
oligomeric enzyme sequentially or independently as already discussed in the Introduction 
section. In the oligomeric enzyme kinetics reaction, the substrate molecules bind to the 
subunits of the oligomeric enzyme in a stepwise manner with different affinity which was 
first proposed by Adair to explain the cooperativity phenomenon observed in the oxygen- 
binding to the hemoglobin at equilibrium^. If an oligomeric enzyme consists of nx number 
of homo or hetero type of subunits, then at the chemiostatic condition the substrate-binding 
scheme of the enzyme can be written as, 

K (0) K (i) K (2) k/™- 1 ' k>) k>t-2) k>t-i) 

E ESi ES 2 s ES n ; ES nT _! ; ES ni 

K^ 1 ' K 2 < 2 ) K 2 ( 3 ) K 2 <"' K 2 (™ +1 ) K^T- 1 ) K 2 ("t) 

(1) 

Here ES n represents the conformational state of the oligomeric enzyme in which n number 
of subunits are occupied by the substrate molecules. K^" -1 ) and are designated as the 
total formation and total dissociation rate constants in the n-th reaction step, respectively. 

The above scheme of substrate binding of an oligomeric enzyme can be viewed as a 
generalization of the kinetics of an enzyme having a single subunit given by 

k/ k_ 2 

E + S^ES^E + P, 

k -i V 

which can be further simplified as 

E^ES. (2) 

K 2 

Here Ki = (ki + k 2 ) and K 2 = (k_i + k_ 2 ), are designated as the total formation and total 
dissociation rate constants of ES, respectively. The pseudo first-order rate constants are 
written as ki = k x [S] and k 2 = k 2 [P] where [S] and [P] are the constant substrate and product 
concentration in the chemiostatic condition. Hence the site-dependent total formation and 
dissociation rate constants in the case of the oligomeric enzyme kinetics are similarly defined 
as Ki (n_1 ) = (k! (n - 1} +k 2 (n - 1) ) and K 2 (n) = (k_x (n) + k_ 2 (n) ) where = ki (n-1) [S] and 

=k' 2 (n - 1) [P]. 

The dynamics of the substrate binding mechanisms are quantified by counting the number 
of occupied sites present in the oligomeric enzyme at a particular instant of time. If at time 
t, 'n' number of occupied sites are present in the oligomeric enzyme (the state ES n ) then 
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a) Kinetic scheme of spatial cooperativity 
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b) Kinetic scheme of temporal cooperativity 
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FIG. 1. Kinetic schemes for (a) spatial and (b) temporal cooperativity with site-dependent binding 
and dissociation rate constants. The numbers in the square boxes denote the number of occupied 
sites. In spatial cooperativity (a), the forward and backward transition probabilities depend only 
on the total formation and dissociation rate constants, respectively. For temporal cooperativity (b), 
the forward transition probability depends on the total formation rate constant and the number 
of unoccupied sites, whereas the backward transition probability depends on the total dissociation 
rate constant and the number of occupied sites. 

at time t + dt, the number of occupied sites may be increased or decreased by one unit 
due to the occurrence of a formation or a dissociation reaction. During the time evolution, 
the number of occupied sites is a fluctuating quantity. Therefore, the system performs a 
one-dimensional random walk along the finite number of states where state-n of the system 
is equivalent to the conformational state ES n , as shown in figured! 

The kinetic scheme of the spatially cooperative system is given in figure QJa). As the 
subunits get occupied sequentially starting from subunit-1, so when we say that the n-th 
subunit is occupied it automatically means that n number of sites are occupied in total. Here 
the forward and the backward transition probabilities depend only on the total formation and 
dissociation rate constants and , respectively which are generally site-dependent. 
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This is so because after the filling of one subunit, there is no other choice for the next 
substrate molecule but to fill up the next adjacent subunit and as this is true for all the 
subunits, there is no combinatorial term in the transition probability. 

The kinetic scheme for the temporal cooperativity is shown in figure QJb). Here the 
substrate molecules can bind independently with any one of the nx number of subunits. 
The state-n of the system represents n-number of occupied sites of the enzyme. In this 
mechanism, the forward transition probability of the n-th state at time t is given by the 

(n) 

product of the total formation rate constant with the number of distinct combinations 
of unoccupied sites present at that time. Similarly, the backward transition probability of 
the same state is the product of the total dissociation rate constant, and the number 
of distinct combinations of occupied sites present at time t (see figure [H(b)). Here these 
rate constants are taken to be site-dependent in general. If all the rate constants are site- 
independent, then the system will be non-cooperative. The main difference of the sequential 
binding scheme from the independent one is as follows: for the sequential scheme, the system 
will show spatial cooperativity in substrate binding even when the formation and dissociation 
rate constants are not site- dependent. 



B. Master equations 

For the time-dependent description of the spatial cooperativity, we have constructed the 
corresponding master equation for this cooperativity mechanism as 

dPsP ^ t] = K^^n - 1; t) + 4 n+1) P sp (n + 1; t) - (K< n) + 4 n) )P sp (n; t), (3) 

with K^ 1 ^ = Kg = K^ nT ^ = Kg = to match the boundary terms. Here, P sp (n, t) is 
the probability of having n number of occupied sites at time t. We have given an analytical 
expression for the solution of the master equation by setting d p sp( n >t) = 0. The steady state 
distribution of the spatial cooperativity is given by 



Psp(n) = " J 7L-i^ . (4) 



where = Ja+i) = k a+i)^a+i) with j = 0, 1, (nr — 1). Here the steady state is actu- 
ally a nonequilibrium steady state (NESS) as already discussed. If X.® = XVj, the NESS 



probability distribution becomes 

P sp( n ) = I _ x(n T +l) ' ( 5 ) 

which is a geometric distribution. The average population of the occupied sites at the NESS 
for X^) = XVj is given by 

/ \ \T- t,«./ x X(l - (n T + l)X nT +n T X nT+1 ) 

(n) = > nP ss fn) = — — ' — -. f6) 

V 1 / Z^ 11L spK ll J (1 - X)(l - X n T+!) K J 



For temporal cooperativity, the master equation is written as 

d Ptem ^ (ll,t) = Kf- 1} (n T - n + l)P temp (n - 1; t) + Kf +1) (n + l)P temp (n + 1; t) 

- K[ n) (n T - n)P tcmp (n; t) - l4 n) nP tcmp (n; t), (7) 



again with K[ = = K^ nT ^ = K^ T+1 ' = 0. Solving this master equation at the NESS, 
we can obtain the probability distribution as 



lj=0 

n " 

K,»M = , (8) 

n I 



c (j) 

occupied sites at the NESS is simply expressed as 



where Xy> = — ottj as aheady defined with j = 0, 1, (np — 1). The average number of 




En 11 

(n> = . (!)) 




Now positive cooperativity in this scenario means a higher affinity of a second sub- 
strate molecule to attach with the oligomeric enzyme compared to that of the first substrate 
molecule which is already bound and so on. Therefore, in this case, the successive bind- 
ing affinity of the substrate molecule increases. So naturally here we take the binding rate 
constants, as follows^^ 

kS 0) < k« ... < kS n) < kS n+1) < .... < k^- 1 ). (io) 



s 



Here the site-dependent overall association rate constant is denned as = H-k^ 
and similarly, the overall site-dependent dissociation rate constant is written as = 
k^{ +k ( i2- We take the rate constants k_j, kg and k^ to be site- independent. This is due 
to the fact that to get a cooperative behavior for the independent binding case, it is not 
necessary to take all the rate constants of the reaction system to be site-dependent that will 
also make the results obtained hard to analyze. Then the site-dependent quantities X.® for 
positive cooperativity maintain the relation: 

x(°) < x«. ... < x( n ) < x( n+1 ) < .... < x^- 1 ). (11) 

Similarly, negative cooperativity arises as a second substrate molecule binds to the 
oligomeric enzyme with a lower affinity than that of the first substrate molecule. Therefore, 
the substrate binding reaction rate constants for different sites obey the inequalities 

kS 0) > k«... > kl n) > kS n+1) > ... > (12) 

Then taking the rate constants k^ , k^ and k^j as site- independent constants, we have 

X<°> > X«. ... > X< n > > X< n+1 ) > .... > X^" 1 ). (13) 

If all the association and dissociation rate constants are site- independent, then the enzyme 
becomes non-cooperative. Thus the steady state distribution, Eq.© reduces to 

p!!( '^(':;)(t^r < i4) 

where X = ^ . By inserting the value of X, the above equation can be written as a binomial 
distribution given by 

^-(:)(^)"(i^r-^- 

This is expected, as in the absence of any cooperativity, the distribution of the occupied 
sites must follow a binomial distribution. So for a system with no cooperativity, the average 
number of occupied sites at the NESS is 



and the average number of unoccupied sites is 

<n T - n) = n T (j^^r) • (17) 

We mention that, in addition to the overall association and dissociation rate constants being 
site- independent, if the rate constant k2 is also negligibly small, then the enzyme kinetics 
becomes simply the Michaelis-Menten type. If (j = 0, . . . , (up — 1)) is taken to be much 
less than the other rate constants, then we have 

X<» = JL, (18) 
where = — 1 2 can be described as the stepwise Michaelis-Menten constant. 

C. Entropy production rates 

The system entropy is defined in terms of the Shannon entropy as 

S sys (t) = -k B ^P(n,t)lnP(n,t), (19) 
ii 

where P(n, t) is the probability of having n number of occupied states at time t with P(n, t) = 
P sp (n, t) or P(n, t) = P te mp( n , t). Here we set the Boltzmann constant, k B = 1. Using the 
master equation, one can get the system entropy production rate^r— as 

Ssys(t) = - ^[w^(n - z/ At |n)P(n - t) - w_ M (n|n - ^)P(n, t)] 

x In^f^i. (20) 
P(n,t) 1 ' 

Here the state of the system can change by any one of the four reactions, denoted with index 

fi, via which the substrate and the product molecules can bind with the enzyme sites and 

k$ n > ._. __ . kL n ? 



detach. They are given as: (1) (ES n + S) 4 (ES n+ i) (// = 1), (2)(ES„) -4 (ES n _i + S) 
{H = -1), (3)(ES n ) ^ (ES n _! + P) (/i = -2) and (4)(ES n + P) % (ES n+ i) (/i = 2). Here 
Ufj, is designated as the stoichiometric coefficient of the /i-th reaction with rate constant k^ 
where = 1 with fi > and — = 1 with // < 0. The transition probabilities are defined 
as follows 

w> - z/» = k(f-^)(n T - (n - !/„)), m > 
10 



and 

w> - z/» = kjr^(n - Vy ),n< 0. (21) 

We have assumed ideal reservoir (surroundings) with no inherent entropy production ex- 
cept through the boundaries of the system. The system entropy production rate(epr) can 
be split as^ 5 - 

S sys (t) = S tot (t)-S m (t). (22) 

Here the first term in the r.h.s. of equation( l22l gives the total entropy production rate and 
the second term denotes the medium entropy production rate due to the entropy flux into 
the surroundings. Therefore the total and medium entropy production rates are defined as 

Stot(t) = - £[w M (n - z/ M |n)P(n - t) - w_^(n|n - z^)P(n, t)] 



2 



x ln W,(n-^|n)P(n-, t) 
w_ M (n|n - f M )P(n,t) 



and 



Sm(t) = - £[w M (n - 2^|n)P(n - z^, t) - w_ M (n|n - ^)P(n, t)] 



2 

x ln ^4^>L. (24 ) 
w_ M (n|n - Vp) 

At steady state, S sys = (whether equilibrium or NESS). An NESS is characterized by a 
non-zero total epr given by 

Slot ESS) = £ [wi(n - l|n)P(n - 1) - w_x(n|n - l)P(n)] 



X ln Vw_i(n|n - l)P(n) x w 2 (n - l|n)P(n - 1) ) ' ^ 
This equation is derived using the circular balance condition 57 

Wi(n — l|n)P(n — 1) — w_i(n|n — l)P(n) = 

w_ 2 (n|n - l)P(n) - w 2 (n - l|n)P(n - 1). (26) 
Now here we consider two limiting situations. 

(i) It is clear that if we do not consider the presence of the product species then there 
will be just two sets of rate constants, Iq and k_l- Then at the steady state, the balance 
condition that holds is obviously the detailed balance which gives 

wi(n - l|n)P(n - 1) - w_i(n|n - l)P(n) = 0. 
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Then from Eg. (125]) . we have S to t = at the steady state which is reduced now to an 

equilibrium. Now from Eq. (ll8[) . the quantity X^ under this condition becomes X^ = kl (j+ ^ 

k -i 

k 'a) 

where -^j+rj are the stepwise equilibrium (binding) constants. So in this limit, theoretically 
there is no difference between our model and a protein-ligand binding model which generally 
does not consider the product formation. 

(ii) Another interesting point is that if we consider the case where ,k_] >> k_2 
(with kj, 1 ^ being already considered negligible) then this corresponds to the pre-equilibrium 
limit or simply the equilibrium limit. The assumption is valid when fast reversible reactions 
precede slower reactions in a reaction network. Now in this situation also, the quantity 
Xu) is defined in terms of the stepwise equilibrium (binding) constants. In this context, 
we mention that the original derivation of the enzyme catalysis reaction by Michaelis and 
Menten involved the pre-equilibrium assumption with the equilibrium dissociation constant 
parameter. The more general derivation by Briggs and Haldane used the steady state ap- 
proximation and their expression contained the actual Michaelis- Menten constant. In our 
case also we see the same features in the quantity, X^ which is the parameter of our model 
study. In the general nonequilibrium case, X^ is related to the stepwise Michaelis-Menten 
constant, (see Eq.( !T8|) with k^ considered negligible) whereas in the absence of product 

species leading to equilibrium or under the pre-equilibrium assumption, X^ is related to the 

k 'a) 

stepwise equilibrium (binding) constant, ^$+t>- 

III. NUMERICAL SIMULATION OF ENTROPY PRODUCTION 

In this section, we have calculated the medium, system and the total entropy produc- 
tion for the spatial and temporal cooperative systems. For a given initial condition, the 
oligomeric enzyme system reaches NESS at a particular time which depends on the chemio- 
static condition, i. e, the value of the constant substrate concentration. The initial condition 
is taken as the fully unbound state of the enzyme with all the subunits being vacant i.e., 
P(n,t = 0) = 5 n ^. This condition leads to zero system entropy at t=0. For the time- 
dependent system entropy production calculation in general, one needs the time- dependent 
solution of the master equation, P(n, t). But here the final time in the calculation of the 
entropy production over the time interval (starting at t = with the specified initial con- 
dition above) is taken such that by then the system reaches the NESS and hence steady 
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state solutions are all we need to get the system entropy production over the length of the 
trajectory. The total entropy production for a single trajectory is calculated over the time 
interval where the determination of the medium entropy production requires the detailed 
information of the path and not just the initial and final points. We run the simulations in 
all the cases up to a fixed point of time taken to be the same for all the binding mechanisms. 
As the steady state is an NESS (and not an equilibrium), total and medium entropy pro- 
duction increase linearly with time and hence if the final point of time is not the same for 
all the cases, one can not compare the various entropy production values for the different 
cooperative systems. 

A. Implementation of the scheme of single trajectory stochastic simulation 

Along a single stochastic trajectory the system entropy production can be defined as^l 
S(t) = — lnp(n, t), where p(n, t) is the solution of the stochastic master equation for a 
given initial condition, p(n ,t ), taken along the specific trajectory n(t). Note that, the 
single trajectory entropy is denoted by (bold) S whereas the trajectory- average entropy 
production (equivalent to ensemble average) is denoted by S. Now at the microscopic level, 
the number of occupied sites of the oligomeric enzyme becomes a fluctuating quantity due 
to the random occurrence of the different reaction events within the random time interval. 
This develops the concept of different trajectories. Here the state of the system can change 
by any one of the four reactions (denoted with index /z) as discussed in Sec.IIC. 

A stochastic trajectory, n(t) starting at the state no, jumping at times tj from the state 
nj_! to the state nj and finally ending up at ni with t = t\ is defined as, 



Here nj = nj_i + where ujf is the stoichiometric coefficient of the /z-th reaction along 
a trajectory and tj = tj_i + 7] where 7] is the time interval between two successive jumps. 
During the jump from the (nj — 1) state to the nj state, any one of the four reactions will 
occur (see Eq. (TjQ) and Eq. (T2])). The rate constant of the reaction /i is denoted as k^. 
The time interval 7] between the two jumps is a random variable following the exponential 
distribution 49,50 



V (!) U W 

n(t) = (n ,t ) -»■ (ni,ti) -> 




(27) 



with a = 



±2 

At=±l 



p(rj) = a exp(— aTj) (28) 
w(nj — 1; ujj. Here w(nj_i; ujr) denotes the forward transition probability 
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from the state (iij — 1) to the state rij through a reaction channel /i with the stoichiometric 



coefficient i/jP . 

Now a time reversed trajectory can be defined as, 

n R (t) = (n b ti) H (n H) tn) ) ... -> (n j; tj) ~% (n^x, tj_i)... (n x ,ti) ^ ' (n ,t ) 

(29) 

This time reversed trajectory is generated due to the occurrence of a reaction channel whose 
state changing vector —ujP is exactly opposite to the state changing vector ujP of the forward 
reaction channel. 

The time-dependent total entropy production, AS to t along a trajectory can be split into 
a system part, AS sys and a medium contribution, AS m . Hence the change of total entropy 
along a trajectory can be written as^- 

AS tot = AS m + AS sys (30) 

where 

AS sys = In^l (31) 
p(n,t) 

and 

ASm .y>?^jfl ,32, 

j w(nj 

Here w(n,-_i; ujP) denotes the forward transition probability as already defined. Similarly, 
w(nj; —vjP) denotes the backward transition probability from the state nj to the (m — 1) 
state through a reaction channel \i with the exactly opposite stoichiometric coefficient —v$. 



B. Cooperative kinetics 

To simulate the spatial cooperativity associated with the sequential binding, we have 
taken the site-independent reaction rate constants as k' : = 0.015 /iM~ 1 s _1 and k_ x = 7.0, 
k_2 = 2.0, k2 = 0.001, all in s _1 . The substrate concentration is taken in fiM unit. The total 
number of subunits present in the oligomeric enzyme is taken as = 3. We have calculated 
the various entropy productions using the stochastic simulation for single trajectories over a 
time interval starting from the initial condition to a final time as mentioned above. We have 
taken 2 x 10 5 trajectories to get the ensemble average of the entropy production values. We 
have calculated the average binding number, (n) for this case from Eq.([6]) and the net product 
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FIG. 2. (a) (n) and (b) v net for the spatial cooperative (sequential) binding as a function of 
substrate concentration, [S] (in fjM. unit) at the NESS, (c) and (d) exhibit the corresponding AS to t 
and AS m variations with [S]. The entropy productions are calculated over a time interval that 
starts with the given initial condition (see text) and ends with the system at the NESS. 

formation rate by using the formula, v net = k_ 2 (n) — k 2 (nx — n), at the final time where the 
system resides at the NESS. We have plotted these quantities as a function of the substrate 
concentration in figure [2£a) and (b). It is clear from the plots that both the quantities 
grow with a sigmoidal shape as a function of substrate concentration indicating positive 
cooperativity in substrate binding. According to Eq.Q, this is due to the higher power 
(> I) dependence of (n) on the factor X which is proportional to the substrate concentration. 
As the rate constants are taken as site-independent, the positive cooperativity generated in 
the system is inherent in the binding mechanism. Now we have plotted AS to t and AS m , both 
being ensemble averages taken over the 2 x 10 5 realizations of the trajectories, in figure [2](c) 
and (d), respectively, against the substrate concentration. Interestingly, we find the nature 
of both the curves to be sigmoidal. 

Next we come to the case of independent substrate binding that can give rise to the 
case of temporal cooperativity with site-dependent reaction rate constants. To simulate 
the entropy production for the positive cooperative system, we take the rate constants of 
successive substrate binding steps as (see 

Eq.<m)- k! 1] = f (1) ki 0) and k< 2) = f^, where 
kj 0) = k} 0) [S] with k} 0) = 0.015 /dVHs- 1 . The set {k} 0) , k_i, k_ 2 , k 2 } is called the starting 
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FIG. 3. (a) (n) and (b) v not for the temporally cooperative (independent) binding with positive 
cooperativity against substrate concentration, [S] (in fjM. unit) at the NESS, (c) and (d) give the 
corresponding AS tot an d AS m variations with [S]. The entropy productions are calculated over a 
time interval as described in the caption of figure [2j It is evident from the figure that all the curves 
show a significant sigmoidal behavior indicating the positive cooperativity. 



or initial rate constants of the cooperative system. For the simulation, here we take f ^ = 10 
and f( 2 ) = 100, i.e., a 10-fold increase in substrate binding rate constants in each step. The 
other rate constants are site-independent and taken to be the same as in the case of the 
spatial cooperativity. We also calculate the average binding number using Eq.(j9]) and the 
net product formation rate, at the NESS. They are shown in figure [3] (a) and (b) along with 
the total and the medium entropy production in figure ^c) and (d), respectively, all as a 
function of the substrate concentration. It is evident from the figure that all the curves show 
a significant sigmoidal behavior indicating the positive cooperativity. We have also given the 
corresponding quantities in the case of non-cooperativity in the same plot for comparison. 
The non-cooperative case is simulated with site-independent rate constants same as in the 
case of the spatial cooperativity. In this case (n) is determined using Eq.f lT6|) . We see that 
in this case also, the nature of variation of (n), v net , AS m and AS to t with the substrate 
concentration is the same, hyperbolic to be specific. 

Now we come to the last case in this category, i.e., the negative cooperativity. In this case, 
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the rate constants of the substrate binding reaction are taken as (see Eq. (|12p ): k/ ^ = 1.5 

/1M-V 1 , kj 1} = f (1) kj 0) and kf } = f( 2 )kj 0) with the values of the factors being f« = 0.1 

and f*- 2 -* = 0.01, i.e., a 10-fold decrease in substrate binding rate constant in each step. The 

other rate constants are taken as in the previous cases. The value of k x is taken to be 100 

times greater compared to the cases of spatial and positive cooperativity. This is only for 

the demonstration of the negative cooperativity effect on the binding curves of the reaction. 

We have plotted (n) against substrate concentration in figure H](a) for the negative as well as 

' (o) 

the non-cooperative case. Here for the non-cooperative case also we have taken k x = 1.5 
//M _1 s _1 . Both the curves show the hyperbolic nature. The two cases are distinguished 
by plotting -A- versus w which is the Lineweaver-Burk plot. For non-cooperative enzyme, 
this plot gives a straight line whereas the curve for the negative cooperative binding starts 
at a higher value on the y-axis and becomes nonlinear when it comes close to the curve 
of the non-cooperative system at high substrate concentration. This feature is shown in 
figure 11(b). Now we have plotted similar curves for AS to t m figure H(c) and (d). One can 
see the same hyperbolic nature in the plot of AS to t versus substrate concentration (figure 
[4(c) ) for both the cases and the nonlinearity in the plot of A ^ - versus at high substrate 
concentration for the negative cooperativity (figured)). So from the above discussion and 
the plots, we conclude that the familiar indications of the cooperative behavior in substrate 
binding, given in terms of the nature of variation of the average binding number and the 
net velocity of the reaction as a function of the substrate concentration, are all reflected in 
the same manner in the corresponding variation of the total as well as the medium entropy 
production. 



We have also calculated the total entropy production rate, S to t at the NESS using Eq. (|25p 
for all the cases of cooperativity . Here we have taken the same set of rate constants as we 
have already considered to calculate the various entropy productions. The variations of S to t 
with substrate concentration, [S] for different binding schemes are shown in figure [5] It is 
evident from the figure that the features of cooperative binding are also reflected in a similar 
fashion on the variation of S to t with substrate concentration. 

17 




FIG. 4. Plots of (a) (n) vs. [S] (in /iM unit) and (b) -X vs. ij at the NESS, (c) AS tot vs. [S] and 
(d) A g t - versus w for negative cooperative (temporal) as well as non-cooperative binding. The 
entropy productions are calculated over a time interval as described in the caption of figure [2j 




FIG. 5. (a) Plot of Stot against substrate concentration, [S] (in /iM unit) for spatial cooperativity. 
In (b) and (c), the same quantity is plotted for positive and negative cooperative cases, respec- 
tively. The non-cooperative case is also shown for comparison, (d) Plot of vs. jgj, which is a 
Lineweaver-Burk type plot, for negative and non-cooperative cases. 
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FIG. 6. Plot of AS sys against substrate concentration, [S] (in //M unit) for (a) spatial cooperative 
binding, (b) positive (temporal) cooperative binding, (c) negative (temporal) cooperative binding 
and (d) non-cooperative binding. In all the cases, the final time of observation is the same, set as 
such that the system reaches the NESS. 

C. System entropy production and binding characteristics 

The ensemble or trajectory average system entropy production over the time interval can 
be written as 

11 T 

AS sys = Sgf - S^ ial = - P ss (n)lnP ss (n), (33) 

n=0 

where the initial condition (time t = 0) of the fully unbound enzyme gives S^ 13,1 = 
and the final state of the system is an NESS characterized by the distribution P ss (n). We 
have plotted the ensemble average system entropy production, AS sys as a function of the 
substrate concentration in figure [6]for all the cases. In figure[6](a), AS sys is plotted for spatial 
cooperativity and in figure [6](b-d) it is shown for the positive, negative and non-cooperative 
cases, respectively which belong to the class of temporal cooperativity. The first thing 
evident from the plots is that AS sys passes through a global maximum for all the cases and 
in the case of negative cooperativity, there is also a local maximum with the parameters of 
our system. 

We have plotted P ss (n) as a function of the substrate concentration in figure C3a-d) with 
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the steady state ('ss') superscript being dropped for simplicity. Figure [7(a) shows the curves 
for spatial cooperativity. We can see that they all cross almost exactly at the same point, 
[S] ~ 600 fiM giving rise to the maximum in AS sys for spatial cooperativity at this point 
(see figure E(a)). In figure [7(b), we have shown the curves for the positively cooperative 
system. At [S] ~ 60 fiM, the curves cross in a pairwise fashion; curves of P(0) and P(3) 
cross each other at this point as well as curves for P(l) and P(2). This particular substrate 
concentration corresponds to the maximum of AS sys in this case (see figure [6(b)). 

The case of negative cooperativity requires a bit more attention. There is again pairwise 
curve crossing of the two sets of probabilities same as in the case of positive cooperativity at 
the same substrate concentration shown in figure[7(c). This gives rise to the global maximum 
in the curve of AS sys for this type of binding shown in figure [6(c). The local maximum can 
be explained as follows. Unlike the plots in figure [7(a) and (b), the probability curves P(2) 
and P(3) remain at significant values over the substrate range studied and the dominance 
of these two probability curves in figure [7(c) (actually when P(2) and P(3) cross, they are 
close to 0.5 at [S] ~ 1800 fiM) over a large substrate range gives rise to an increase of 
AS sys , albeit slow. Finally we come to the case of non-cooperativity in figure [7(d) where 
again there is the pairwise crossing of the same set of probabilities as in figure [7(b) but 
at [S] ~ 600 /iM that again gives rise to the maximum of AS sys shown in figure [6](d). In 
this case too, there are more than one dominating probability curves before and after the 
pairwise crossing over similar substrate range as in figure [7(c). But the AS sys in this case 
shows a slow but steady decrease with substrate concentration after passing through the 
maximum without any unusual behavior. This may be due to the fact that here at least 
three of the four probabilities are significant (with comparable values) over a large substrate 
range and so they do not cross the value of 0.5 in this range unlike the case in figure [7(c). 
It is clear that arbitrary variation of the rate constants of the system in each binding step 
can make life more complicated and then the maxima in the AS sys curve may or may not 
be associated with the binding probability curve crossings. 

We can justify the curve crossings, whether they all cross or cross pairwise at a particular 
substrate concentration, by inspecting the expressions of the steady state probability dis- 
tributions. We see from the steady state distribution for the spatial cooperativity, Eq.([5]), 
that if one of the probabilities, say P(0) is approximately equal to any other probability, say 
P(3), then X ~ 1 (but obviously not exactly equal to 1) and this automatically leads to the 
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FIG. 7. Plot of the steady state distributions, P ss (n) against [S] (in fj,M unit) for (a) spatial coop- 
erative binding, (b) positive (temporal) cooperative binding, (c) negative (temporal) cooperative 
binding and (d) non-cooperative binding. In the plots, the 'ss' superscript is dropped for simplicity. 
The arrows indicate the curve crossing points. 

near equality of all the probabilities at this value of X. This is true for all the probabilities 
and hence in this case the probabilities can only cross simultaneously at X ~ 1. As the prob- 
abilities are equal at this point which corresponds to [S] ~ 600 fiM, this obviously gives the 
maximum system entropy production in this case. Now we take the steady state distribution 
of the non-cooperative case, Eq. ffT4"j) . It can be easily seen that here only P(0) = P(3) leads 
to the equality P(l) = P(2) at X = 1 giving the maximum of AS sys again at [S] ~ 600 /xM. 
So in the context of the system entropy production the spatial cooperative system shows 
some similarity with the non-cooperative system. With the same set of site-independent 
rate constants, the spatial cooperative system is also associated with larger system entropy 
production compared to that of the non-cooperative case. This is because all the binding 
probabilities become equal for the spatial cooperativity whereas they become equal pairwise 
for non-cooperative binding. 

The cases of positive and negative cooperativity belonging to the class of temporal co- 
operativity are a bit complicated. We have considered a 10-fold increase of the substrate 
binding rate constant for each successive binding in the case of positive cooperativity whereas 
a 10-fold decrease in the corresponding rate constant for each successive binding for negative 
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cooperativity. This symmetry ensures that in both the cases only P(0) = P(3) leads to the 
equality P(l) = P(2) at X« = 1. This can be easily proved from Eq.flH]). But if the rise or 
fall of the value of the substrate binding rate constant in each successive step of binding is 
not by the same factor, then the pairwise equality of the binding probabilities is not possible 
at a given substrate concentration. 



IV. MEASURE OF COOPERATIVITY 

Here we have discussed on the determination of the Hill coefficient from the master 
equation corresponding to the different binding schemes. We have also introduced an index 
of cooperativity in terms of the stochastic system entropy associated with the fully bound 
state of the cooperative and non-cooperative cases. We have analyzed its connection to the 
Hill coefficient using some relevant experimental data which gives a realistic application of 
the proposed scheme of measurement of cooperativity. 



A. Hill coefficient 

In the traditional enzymology, the characterization of cooperativity is carried out by 
measuring the Hill coefficient^. For positive and negative cooperative cases, the Hill co- 
efficient becomes greater than or less than one, respectively, whereas the non-cooperative 
case is characterized with Hill coefficient equal to one. Experimentally it is obtained by 
determining the fractional saturation, 9(= (n)/nx) at various substrate concentrations [S], 
constructing the Hill plot (ln(j-zg) vs. ln[S]) and then finding the slope at the half-saturation 
point, 9 = 0.5 or at a point where the slope deviates maximum from unity. On the other 
hand, Hill coefficient is theoretically defined as the ratio of the variances of the binding 
number of the cooperative and non-cooperative cases at the half-saturation point with the 
non-cooperative binding case following the binomial distribution^^. 

We briefly mention the features of the Hill plot for the model binding schemes studied 
here. The slope of the Hill plot is generally given by^ 

[S](d0/d[S]) 

H = 9(1-9) ■ (34) 
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For temporal cooperativity, the fractional saturation can be written as (see Eq.(j9])) 

q _ En"B n [Sf 



where B n = ( n ^) n^o^S) -1 with B = 1. Then one gets 

(n 2 ) - (n) 2 _ (T t 2 emp 

n t emp - fth m - 2 ' 

n T 0(l - 0) <r bino 

where cx 2 emp and o" 2 ino are the variances of the binding numbers of the temporal and non- 
cooperative cases, respectively. The Hill coefficient, nn is given at the half-saturation point 
as^ nn = — ^— • Similar expressions hold for the spatial cooperative binding. We have 
already mentioned in Sec.IIB that if all the rate constants of the independent binding scheme 
are site- independent, then the binding is non-cooperative with binomial distribution of the 
binding probability. Here we discuss the corresponding scenario for the sequential binding 
(leading to spatial cooperativity) to be non-cooperative in terms of the variance of the 
binding number. In the case of spatial cooperativity, the variance of the binding number, 

°s 2 P is S iven b Y 

2 /1 + Xn 2 n T (n T + l)X n *+i 

% = [T^x) {n) ~ (n) 1 _ x*t+i — ' (37) 

where X and (n) are as given in Eq.(j5]) and Eq.flB}. Now for n-r = 1, this reduces to 

^sp = (! +X )2 = ^ino- (38) 

2 

and then the slope of the Hill plot becomes H sp = = 1 for any substrate concentration. 
So for sequential binding, the cooperative behavior is absent only if the enzyme is monomeric. 



B. Cooperativity index 

Here we introduce an index of cooperativity. First we build up the concept from binding 
probabilities and then demonstrate how this index can indicate the nature of the cooperativ- 
ity. For positive cooperative binding, one expects that full occupancy of the enzyme is more 
probable compared to the case of non-cooperativity. Similarly, for negative cooperativity, 
the full occupancy of the enzyme is less probable. Now, if the probability of an event-n is 
p n , then the associated surprisal is given by — ln(p n ) and more probable the event, the less is 
its surprisal. So the ratio of the surprisals, associated with the probability of the system to 
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remain in a fully occupied state without and with cooperativity at NESS, should be greater 
than 1 for positive cooperativity and less than 1 for negative cooperativity. Hence we define 
the index of cooperativity, denoted by C in terms of the ratio of the surprisals, associated 
with the probability of the system to remain in a fully occupied state without and with 
cooperativity at NESS as 

_ -l n ( P (bino) (nT)) 

° " -In(Q(n T )) (39) 
where the binomial distribution, p( bino ) is the reference corresponding to the non-cooperative 
case and the distribution Q corresponds to the cooperative binding case. The rate constants 
of the reference non-cooperative system (binomial) must be the same as those of the starting 
or initial rate constants of the cooperative system for the comparison to be valid. The relation 
is then independent of the actual value of the (constant) substrate concentration. We point 
out that the surprisal is equivalent to the single trajectory stochastic system entropy2i>£2 
(associated with the fully occupied state). So the index, C is truly an entropic estimate of 
cooperativity at the microscopic level. 

Based on the above argument, next we theoretically analyze the probability to remain 
in a fully occupied state for different cooperative systems and compare that with the non- 
cooperative case to formulate the criteria of cooperativity in terms of C. For spatial co- 
operativity, the ratio of its steady state distribution (Eq.(j5])) and the reference binomial 
distribution (Eq.f JT4|) ) for n = nx is given by 



R. 



"P _ p(bino) | 



n T J 



1 + 



((7) - ( n T )) x + ((7) - (7)) x2 + • • • + (© - ( nT nT -i)) xnT 



(40) 



(1 - X n T+!) 

^From the above expression it is clear that for all values of X, either greater than or less 
than 1, the quantity R sp is greater than 1 indicating positive cooperativity. This will lead 
to the condition C > 1 for the case of spatial cooperativity for any substrate concentration. 
Only in the case of monomeric enzyme, (n^ = 1), the system will be non-cooperative with 
Rsp = C = 1 as already discussed in terms of variances at the end of Sec.IVA. 

In the case of temporal cooperativity, the corresponding ratio, Rt em p is given using Eq.® 
and Eq. (fT4l) at n = n T as 

^ Ptemp( n T) 

rl 



"'"'!> - p(bino)( nT ) 
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x(°)x( 1 )...x( n T- 1 ) 

[l+n T X( )+...+(x'( )X( 1 )...X( n T- 1 ))] 

X"T 



(l+X) n T 

(X(°') n Tf( n T- 1 ) 
[l+Il T X( )+... + (X(°)) n Tf("T- 1 )] 
(X<°)) n T 
[l+n T X(°)+...+(X(°)) n T] 



(41) 



C ( 



with X(°) = X. Now both P^ mp (n T ) and p( bino )(n T ) tend to 1 at large X<°> i.e. large substrate 
concentration. But it is clear that for positive cooperative binding with f > 1, the last term 
in the denominator of Ptemp( n T) dominates the previous terms more readily compared to 
the case of p( bmo )(nr). Hence at a particular substrate concentration, PtI m p( n T) is closer to 
1 compared to P( bino )(nT-) and so Rt em p is greater than 1. For negative cooperativity with 
f < 1, the situation is obviously reverse and Rtemp is less than 1. Therefore, in the light of 
the above discussions and Eq. fl39|) . we write down the condition of cooperativity in terms of 
C as 

> 1, positive cooperativity 

= 1, no cooperativity (42) 
< 1, negative cooperativity. 
This is the same criteria of cooperativity as given in terms of the Hill coefficient. To find out 
the Hill coefficient, i.e., the variances theoretically, it is necessary to know the probability 
distribution of the corresponding positive and negative cooperativity cases, respectively. 
Now our measure of cooperativity, the index C, is also related to the probability distributions; 
but it is defined in terms of the ratio of a specific term of the distributions, namely the 
probability of the fully occupied state. So apparently there is no straightforward connection 
between the Hill coefficient and C. The Hill coefficient is the slope of the binding curve at 
a particular substrate concentration corresponding to the half-saturation point whereas the 
index C is defined independent of the substrate concentration and the characterization of 
cooperativity in terms of C is valid at any substrate concentration. 

We have plotted the quantity, C in figure OH]) for positive cooperative system (independent 
binding) and also for the spatial cooperative binding for different values of nx as a function 
of substrate concentration. For the positive cooperativity case, the substrate binding rate 
constants, increase by a factor of 2 in each step. The value of C grows with substrate 
concentration, starting just above 1.0 and finally saturates. One can see from Eq. (|39l . 
that the limiting value of C (obtained at high substrate concentration) in case of spatial 
cooperativity is nx whereas for temporal cooperativity it is given by f( n T -1 ) where f( nT_1 ) = 
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FIG. 8. Plot of the cooperativity index, C against substrate concentration, [S] (in //M unit) for 
different values of the number of subunits of the oligomeric enzyme, np in the case of (a) positive 
(temporal) cooperativity and (b) spatial cooperativity. 
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These results are discussed in detail in the appendix. Here we specifically mention 



the case of np = 5 for the positive cooperativity where the limiting value of C is f ^ 
2 4 = 16. It is evident from figure ©(a) that this is indeed the case. 



4 4) 



C. Characterization of cooperativity: a case study with stepwise Aspartate 
receptor binding 

Although apparently there is no straightforward connection between the Hill coefficient, 
nn and C, first of all it is clear that the well-known criteria of cooperativity in terms of nn 
is exactly the criteria we have given in terms of the cooperativity index, C in Eq. (l42l) . Both 
the measures are equal to 1 in the absence of cooperativity whereas for cooperative binding, 
the criteria are the same although the actual values of nn and C will be generally different. 
Here we will try to illustrate this point using some experimental data from the work of 
Kolodziej et alM regarding the production of positive, negative as well as non-cooperativity 
by mutations at a serine 68 residue located at the subunit interface in the dimeric aspartate 
receptor of Salmonella typhimurium. Due to unavailability of experimental data of the 
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stepwise Michaelis-Menten constants, K^, we use the stepwise binding constants reported 
in their study in the place of (K^) -1 of the independent binding model with nx = 2 
and j = 0, 1. Now the parameter X® in our study is related to in the general non- 
equilibrium condition and reduces to stepwise equilibrium (binding) constants under the 
conditions already discussed at the end of Sec.IIC. For experimental testability of C at 
NESS, one needs the stepwise Michaelis-Menten constants, • We choose the independent 
binding model as the experimental result reports both positive and negative cooperativity. 
We calculate the fractional saturation 9 as a function of substrate concentration, [S] using 
Eq.(|9]) and find out the Hill coefficient, nn at the half-saturation point (6 = 0.5). Then we 
determine the cooperativity index, C at the substrate concentration where 9 = 0.5. The 
results are given in Table. 1. The Hill coefficients derived by us for different cases tally 
very well with the experimental data 52 . The cooperativity index, C detects the presence 
and absence of cooperativity successfully. Also the extent or degree of positive or negative 
cooperative behavior is equally well characterized by the index, C. This can be seen by 
comparing the values of nn and C for the cases of serine and cysteine showing negative 
cooperativity as well as for threonine and isoleucine showing positive cooperativity. 

TABLE I. The stepwise Aspartate binding constants, K 1 and K 2 (in /iM _1 ) for different amino 
acid residues at position 68 of Aspartate receptor taken from the experimental study of Kolodziej 
et a^. Here we have taken the values of the inverse of the stepwise Michaelis-Menten constants, 
Kjjj in our model to be equal to the binding constants. The values of the Hill coefficient, nn in 
the parentheses are from the experimental work, given for comparison with the values determined 
here. The cooperativity index, C characterizes the cooperative behavior successfully as can be seen 
by comparing it with nn- 



Amino acid 


M 




nH 


C 


serine 


0.7 


0.2 


0.7(0.7) 


0.491 


cysteine 


0.5 


0.2 


0.776(0.8) 


0.598 


threonine 


0.4 


0.9 


1.197(1.2) 


1.519 


isoleucine 


0.4 


2.8 


1.446(1.4) 


2.558 


aspartate 


0.1 


0.1 


1.0(1.0) 


1.0 



The cooperativity index, C is related to the probability of fully bound state of the single 
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enzyme. So another possibility of experimentally determining C, apart from the measure- 
ment of the stepwise Michaelis-Menten constants, will be to detect this fully bound state 
by electrical or optical means in a single molecule experiment and then to fit the resulting 
probability with some model distribution. 

V. CONCLUSION 

We have classified the cooperative substrate binding phenomena of a single oligomeric 
enzyme on the basis of the binding mechanism and the nature of the substrate-bound states 
of the system in a chemiostatic condition. Both the binding mechanisms are modelled in 
terms of master equation. The sequential binding of the substrate molecules leads to spatial 
cooperativity whereas the independent binding scheme leads to temporal cooperativity. We 
have determined the various entropy productions due to the enzyme kinetics over a time 
interval where at the final point of time the system is in a nonequilibrium steady state 
(NESS) that can be arbitrarily far away from equilibrium. We have used kinetic Monte 
Carlo simulation algorithm applied on a single trajectory basis to calculate the entropy 
production. In this context, the interesting finding is that the total as well as the medium 
entropy production show the same diagnostic signatures for detecting the cooperativity as is 
well known in terms of the average binding number or the net velocity of the reaction. More 
specifically, AS to t as well as AS m for positive cooperative kinetics show sigmoidal variation 
as a function of substrate concentration whether the class being spatial or temporal. They 
also show the non-linearity in the inverse plot of Lineweaver-Burk type demonstrating the 
case of negative cooperativity. The signs of cooperative behavior is also reflected in a 
similar fashion on the variation of the total entropy production rate (epr) with substrate 
concentration determined at the NESS for different binding schemes. That the features of 
cooperativity are reflected similarly on the variations of both the total epr at the NESS and 
the total (and medium) entropy production over a time interval up to the NESS is a highly 
interesting fact and gives deep insight on the role of the binding mechanism in governing 
the total entropy production in a general non-equilibrium setup. 

We have thoroughly analyzed the system entropy production for all the cases in terms 
of the steady state binding probability distributions. For a spatial and a non-cooperative 
system, the maximum value of the system entropy production due to the nonequilibrium 
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processes in the reaction appears at the same substrate concentration with the value of the 
entropy production being greater for the spatial cooperativity. We have explained this in 
terms of the different binding probability curve-crossings that helps to understand how the 
binding characteristics affect the entropy production of the system, i. e, the single oligomeric 
enzyme. Similarly, the distinct features of the evolution of system entropy production for 
the positive and negative cooperative binding give valuable insights on its connection to the 
binding mechanism. 

We have introduced an index of cooperativity, C defined as the ratio of the surprisal 
or equivalently, the stochastic system entropy associated with the fully bound state of the 
cooperative and non-cooperative cases. The criteria of cooperativity in terms of C is identical 
to that of the Hill coefficient. We have analyzed its connection to the Hill coefficient using 
some relevant experimental data. This index is truly an entropic estimate of cooperativity 
and gives a microscopic insight on the cooperative binding of substrate on a single oligomeric 
enzyme instead of realising cooperativity in terms of macroscopic reaction rate. 
Acknowledgement : K.B. acknowledges the Council of Scientific and Industrial Research 
(C.S.I.R.), India for the partial financial support as a Senior Research Fellow. 



Appendix: Estimate of the limiting value of the cooperativity index, C for 
various cooperative binding 

Here the limiting value of the cooperativity index, C for the spatial and temporal coop- 
erative binding are determined at high substrate concentration. The limiting value of the 
cooperativity index, C for the spatial cooperativity is calculated from Eq. (l39|) by using the 
steady state probability distribution function of spatial cooperativity, Psp(n) (Eq.([5])) and 
that of no cooperativity, Pbj no (n) (Eq.(fT4l)) at n = up. The expression of C then becomes 

C = x) . (A.l) 



_i X n T(l-X) 
111 [l_ X ( n T + l) 

At high substrate concentration, with X >> 1, the above equation can be written as 

n T ln(l + i) 

c = ^4- (A - 2) 



Now expanding the log terms in the Eq. ( ]A.2j) and neglecting the higher order terms, we 
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finally obtain 

C = n T . (A.3) 

Therefore, the limiting value of C in the case of spatial cooperativity, obtained at high 
substrate concentration, is equal to the total number of sub-units of the oligomeric enzyme. 

In a similar fashion, the limiting value of C can be calculated for the temporal coopera- 
tivity from Eq.f l39p by using the steady state probability distribution function of temporal 
cooperativity, Ptg mp (n) (Eq.flH])) and that of no cooperativity, P^ mo (n) (Eq.flHj)) at n = nT- 
At this value the distribution, P^ mp (n) can be written as 

Ptemp(n T ) - [x + ^(0) + + ( X (0)X( 1 )...X( n T-l))] ' (A ' 4) 

Here X (j) « fQ)x (0) with j = 0, (n T - 1). This follows from the definition of X (j) (see 
Eq.flH]) and Eq.(jH)) with the small value of k_ 2 taken in this study. Now at high substrate 
concentration with X^ ) >> 1, the above equation can be written as 

PtcmpK) = F— y (A.5) 

I 1 ~T x ( ) f (n T -l)J 

Now, by using the value of P^ (up) and Pb ino (nr) into the Eq. (l39l) at high substrate 
concentration, we obtain 

C = , , . (A.6) 



ln[l + ^ 



X (0)f(n T -l)J 

For the comparative study of the temporal and non-cooperative cases, the starting value, 
X^ ) is taken equal to X. Then expanding the log terms in the above equation and neglecting 
the higher order terms, we finally obtain the limiting value of C for temporal cooperativity 

as 

i(n T -l) 

C = f (nT " x) = (A.7) 
k i 

Here we mention that for the negative cooperative binding, f( nT_1 ) can be much less than 1 
in general. But here we consider the case X^f^ nT_1 ^ >> 1. 
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